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Abstract: Nitrogen was the main limiting nutrient of net primary pro¬ 
duction in the southeastern Keerqin Sandy Lands, Northeast China. Spe¬ 
cies richness declined and biomass increased after five consecutive years 
of nitrogen fertilization of these sandy grasslands (2004—2008). After 
fertilization had been stopped for three years (2009—2011), we surveyed 
vegetation on previously fertilized plots to quantify changes in commu¬ 
nity composition. Respect species richness showed an increasing trend 
over time since the cessation of fertilization. Respect vegetation height 
and coverage showed decreasing trends over time since the cessation of 
fertilization. Species composition changed after fertilization ceased, the 
dominant species shifting from Cannabis sativa, Phragmites communis 
and Chenopodium acuminatum in 2008 to Cannabis sativa, Phragmites 
communis and Artemisia scoparia in 2011. Dominance of dominant 
species declined from 66.2% in 2008 to 57.5% in 2011. The importance 
value of annual plants in the earlier nitrogen addition plots was higher 
than in control plots, but the differences were not significant in 2011. The 
importance value of perennial plants differed significantly between 
treatments from 2009 to 2011. The reversion rate not only differed be¬ 
tween community characteristics, but also between functional groups in 
the same community characteristic. Although the residual effect of nitro¬ 
gen addition on vegetation was still observed three years after fertiliza¬ 
tion ceased, the vegetation showed signs of recovery. 

Keywords: nitrogen fertilization, recovery, Keerqin Sandy Lands, vege¬ 
tation 


Project funding: The study was supported by the Chinese Academy of 
Sciences, Projects No. XDA05050401 and XDA050570100, and by the 
National Natural Science Foundation of China, Project No. 30800143. 

The online version is available at http://www.springerlink.com 

Sheng-nan Shi u ' 3 , Zhan-yuan Yu u (EiJ), Qiong Zhao 1,3 
1 State Key Laboratory of Forest and Soil Ecology; Institute of Applied 
Ecology, Chinese Academy of Sciences, Shenyang 110016; 2 University 
of Chinese Academy of Sciences, Beijing 100039; ’ Daqinggou Eco¬ 
logical Station, Institute of Applied Ecology, Chinese Academy of Sci¬ 
ences, Shenyang 110016. E-mail: fafuswcl@163.com 

Corresponding editor: Hu Yanbo 


Introduction 

Nitrogen is recognized as the limiting nutrient of net primary 
production in many terrestrial ecosystems (Vitousek et al. 1991; 
LcBauer et al. 2008). Therefore, nitrogen fertilization has been 
widely used to increase plant growth and productivity in terres¬ 
trial N-limited ecosystems (Frink et al. 1999). Increased nitrogen 
input increases soil N availability, which leads to changes in 
species composition and reduces species diversity (Wedin et al. 
1996; Bobbink et al. 1998). 

Biodiversity plays an important role in ecosystem functioning, 
sustainability and stability (Tilman et al. 1994; Tilman et al. 
1996). The effect of nitrogen addition/deposition on species 
composition and diversity has been demonstrated in numerous 
studies (Willems et al. 1993; Stevens et al. 2004; Lu et al. 2010; 
Dupre et al. 2010). In comparison, studies of the resilience of 
vegetation species composition after fertilization was ceased or 
reduced are fewer (Boeye et al. 1997; Giisewell et al. 2002; 
Clark et al. 2010; Edmondson et al. 2013). Resilience is defined 
as the ability of ecosystems to recover after termination of fertil¬ 
izer addition in our study (Leps et al. 1982). To date, the reversi¬ 
bility of N-induced shifts in vegetation is largely unknown 
(Clark et al. 2009). The studies at the Schynige Platte showed 
that short-term and small-scale fertilization perturbation has 
long-lasting effect on subalpine acid grassland, and the ecosys¬ 
tem can recover to its original state with time (Spiegelberger et al. 
2006). Similarly, Boxman et al. (1998) reported that nitrogen 
enrichment was a fully reversible process. In contrast, Semelova 
et al. (2008) reported no signs of plant recovery several years 
after fertilization was terminated. Research conducted in a large 
grassland suggested that recovery of species richness was rela¬ 
tively rapid compared with most other studies (Hegg et al. 1992; 
Walker et al. 2004; Josef et al. 2009). In alluvial grassland a du¬ 
ration of 16 years was adequate to restore species richness, but 
not plant species composition or N/P ratios in plant biomass 
(Hrevusova et al. 2009). O'Sullivan et al. (2011) concluded that 
the recovery of pre-fertilization plant species composition re¬ 
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quired more time than that recovery of soil chemical characters 
in limestone and acidic grasslands. This conclusion applied as 
well to boreal forest ecosystems (Strengbom et al. 2011). These 
studies suggested that the resilience of ecosystems differs not 
only between ecosystems, but also between processes within the 
same ecosystem (Lavorel et al. 1999). The studies were mainly 
conducted in Europe and America where the rate of nitrogen 
deposition had been reduced or stabilized (NADP 2000), on farm 
fields that were previously fertilized on an experimental basis 
(Hejcman et al. 2007) and in artificial roof experiments (Boxman 
et al. 1998). The fertilization/deposition histories of these sites 
differed and this might explain the sometimes contradictory re¬ 
sults. In addition, ecosystem recovery might be influenced by 
factors other than fertilization such as hysteresis, soil seed banks, 
nearby seed sources and other factors. Thus, the studies investi¬ 
gating the vegetation recovery need to be carried out on a global 
scale and over long time periods. 

Keerqin Sandy Lands are the largest of China's four major 
sandy lands. A previous fertilization experiment reported by Yu 
et al. (2006) demonstrated that the limiting factor of the 
semi-arid sandy grassland was nitrogen. The present study was 
conducted on this experimental fertilization site from 2009 to 
2011 after anthropogenic nitrogen inputs ceased in 2008. The 
aim of our research was to test the reversion rate and pattern of 
semi-arid sandy grassland in tenns of plant diversity and com¬ 
munity composition after cessation of fertilization. 

Materials and methods 

Experimental site 

The experimental site was located at Daqinggou Ecological Sta¬ 
tion (latitude 42°58' N, longitude 122°21' E, 260 m.a.s.l.) of the 
Institute of Applied Ecology, Chinese Academy of Sciences, in 
the southeastern Keerqin Sandy Lands of northeast China. The 
area has a typical continental monsoon climate. The mean annual 
temperature is 6.4 °C with minimum and maximum mean 
monthly temperatures of -12.5 °C and 23.8 °C, respectively. 
Mean annual precipitation is 450 mm, with more than 60% from 
June to August. The soil is classified as Typic Ustipsamment with a 
texture of 90.9% sand, 5.0% silt, and 4.1% clay (Zeng et al. 2009). 
The growing season lasts from late April to early September. 

The site was used as cropland in 1997 and then naturally re¬ 
stored to grassland in 2000. Afterwards, the herbaceous vegeta¬ 
tion recovered rapidly to coverage of 80%. The dominant plant 
species were Pennisetum flaccidum, Chenopodium acuminatum, 
and Cleistogenes chinensis. 

Experimental design 

In 2004, a nitrogen addition experiment was conducted on flat 
sandy grassland that was fenced to exclude livestock grazing 
from 2003. We sampled six replicate plots on two treatments 
(nitrogen addition and control). The area of each plot was 16 m 2 
(4 m x 4 m), and there was a 2 -m buffer strip between any two 
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adjacent plots. From 2004 to 2008, N (NH 4 NO 3 ) was added at 
the rate of 20 g N-m' 2 -a''. From 2009 to 2011, fertilizer addition 
was stopped to document vegetation recovery. 

Measurement of plant responses 

Annual vegetation surveys were carried out in mid-August, the 
time of peak biomass, in the period 2008-2011. In each 16 m 2 
plot, four randomly located 1 m x 1 m quadrats were sampled to 
investigate plant abundance and cover. Plant abundance was 
visually estimated according to the Braun-Blanguet scale. The 
cover of each species and the total cover of each plot were esti¬ 
mated directly in percentages. The height of vegetation in each 
plot was estimated using a tape measure. 

Data processing and analysis 

Data processing 

We used plant species richness, the Shannon-Wiener index (H) 
and the Simpson index ( D ) to describe species diversity. Plant 
species richness was the number of different species per square 
meter. The Shannon-Wiener index was calculated as: 

H = —^ Pi InPi, 

The Simpson index was calculated as: 

D = l -IX 

Where P, is the relative importance values of species i, the spe¬ 
cies relative importance value is the mean of relative cover, rela¬ 
tive frequency and relative density. Because climatic variation 
and other factors also affect species diversity, when we analyzed 
differences between years we calculated relative species richness 
by dividing the plant species richness in a given year by the 
mean value in the control plots for that year (Clark 2008). The 
relative Shannon-Wiener index, relative Simpson index, the rela¬ 
tive height and coverage were calculated in the same way. 

Data analysis 

Species richness and species diversity in fonner nitrogen addi¬ 
tion treatment plots and control plots were analyzed using Stu¬ 
dent's t-test. The changes of ecosystem composition after cessa¬ 
tion of fertilization between years were tested using one-way 
ANOVA. LSD method was used for multiple comparisons if the 
data were homogeneous, otherwise Tamhane analysis was used. 
All analyses were perfonned using SPSS 13.0 and the accepted 
significance level was a = 0.05. 

Results 

Plant diversity 

Species richness, diversity and evenness on fonnerly fertilized 
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plots were significantly lower than on control plots during all 
years 2009-2011 (p <0.05; Table 1). Relative species richness 
showed an increasing trend over time since cessation of fertiliz¬ 
ing. Relative species richness increased by 40% (from 0.55 in 
2008 to 0.77 in 2011) and was significantly greater in 2011 than 


in 2010 (p <0.05; Fig. lc). Relative diversity was from 0.75 in 
2008 to 0.89 in 2011. Relative evenness was from 0.87 in 2008 
to 0.94 in 2011. The differences in relative diversity and relative 
evenness between years were not significant (p >0.05). 





Year 


Fig. 1: Changes of relative diversity, relative evenness and relative species richness from 2008 to 2011 (mean±SE, n =6). Different small letters repre¬ 
sented the significant difference of mean values among years (p <0.05). The same below. 


Table 1: Responses in 2009—2011 of diversity, evenness and species 
richness to cessation of nitrogen addition in 2008. 


Year 

Treatment 

Diversity 

Evenness 

Richness 

(m- 1 ) 

2009 

Control 

2.44 (0.04) a 

0.89 (0.01) a 

9.63 (0.50) a 

N 

1.87 (0.13) b 

0.80 (0.03) b 

5.88 (0.59) b 

2010 

Control 

2.36 (0.05) a 

0.89 (0.00) a 

9.08 (0.33) a 

N 

1.86 (0.06) b 

0.79 (0.01) b 

5.83 (0.24) b 

2011 

Control 

2.25 (0.06) a 

0.88 (0.01) a 

8.46 (0.64) a 

N 

1.99 (0.05) b 

0.83 (0.01) b 

6.54 (0.25) b 


Values in parentheses are SE (n =6). Different small letters within a column 
indicate the significant difference of mean values (p <0.05). The same bellow. 

Species composition 

The species composition on formerly fertilized plots differed 
from that on control plots in all years after cessation of nitrogen 
fertilization. The dominant species shifted from Cannabis sativa, 
Phragmites communis and Chenopodium acuminatum in 2008 to 
Cannabis sativa, Phragmites communis and Artemisia scoparia 
in 2011. The dominant species in 2010 were the same as in 2011. 
The dominance of dominant species declined from 66.2% in 
2008 to 57.5% in 2011. Although species composition differed 
by treatment in 2011, C. sativa and P communis were dominant 
on formerly fertilized plots and on control plots. 

Functional composition of species 

The importance value of annual plants on formerly fertilized 


plots was higher than on control plots, the differences were sig¬ 
nificant in 2009 and 2010 (p <0.05 both years) but not in 2011 (p 
>0.05; Table 2). The yearly variation in biennial plants between 
treatments was not significant (p >0.05; Table 2). The importance 
value of perennial plants was significantly lower on fonnerly 
fertilized plots than on control plots in every year after nitrogen 
addition ceased (p <0.05; Table 2). 


Table 2: Responses of functional groups to nitrogen addition ceased in 
2009-2011. 


Year 

Treatment 


Importance Value 

Annual 

Biennial 

Perennial 

2009 

Control 

33.61 (4.05) a 

7.40 (1.56) a 

58.99 (4.64) a 

N 

57.45 (3.79) b 

9.71 (5.14) a 

32.84 (6.62) b 

2010 

Control 

N 

30.78 (2.02) a 

50.84 (5.50) b 

16.33 (2.14) a 

23.92 (6.18) a 

52.89 (2.51) a 

25.24 (2.90) b 

2011 

Control 

39.91 (3.36) a 

9.37 (2.74) a 

50.72 (1.31) a 

N 

48.32 (5.74) a 

14.89 (1.97) a 

36.79 (4.68) b 


Plant height and coverage 

Vegetation height and coverage on formerly fertilized plots were 
greater than on control plots in 2009-2011 (Table 3). The differ¬ 
ences in vegetation height between treatments were significant in 
2009 and 2011 (p <0.05) but not in 2010 (p >0.05; Table 3). In 
contrast, the differences in coverage between treatments were 
significant in 2009 and 2010 (p <0.05) but not in 2011 (p >0.05; 
Table 3). The relative height values in 2010 and 2011 were sig¬ 
nificantly lower than in 2009 (Fig. 2a). Relative height declined 
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by 44% in 2010 and 39% in 2011 compared to 2009. Relative 
coverage declined after fertilization ceased, and coverage was 
significantly greater in 2009 than in 2011 (p <0.05; Fig. 2b). The 
responses of coverage to cessation of nitrogen inputs varied 
among dominant species. The coverage of C. sativa declined 
from 36% in 2009 to 30% in 2011. However, the coverage of P. 
communis increased from 15% in 2009 to 29% in 2011. A. sco- 
paria became a dominant species in 2010 and 2011 when its 
coverage reached 12% and 7%, respectively. 


Table 3: Responses during 2009-2011 of vegetation height and coverage 
to cessation of nitrogen addition in 2008 


Year 

Treatment 

Height (cm) 

Coverage (%) 


Control 

33.33 (2.11) a 

78.33 (2.79) a 

2009 


N 

83.33 (4.94) b 

91.67 (2.11) b 


Control 

65.00 (9.83) a 

85.83 (1.54) a 

2010 


N 

90.83 (12.00)a 

98.33 (1.67) b 


Control 

65.83 (7.35) a 

90.00 (3.42) a 

2011 


N 

100.00 (9.57) b 

98.33 (1.67) a 



Year 

Fig. 2: Changes of relative vegetation height and coverage by year in 
2009—2011 (meaniSE, n =6). 

Discussion 

Species richness, diversity and evenness were dissimilar on for¬ 
merly fertilized plots and control plots. These results showed that 
the residual effect of fertilization was still observed three years 
after fertilization ceased. Inorganic nitrogen concentration in 
soils was significant higher on formerly fertilized plots than on 
control plots in 2009 (Li et al. 2012), but was similar between 
treatments in 2012. We don’t know whether or not the differ¬ 
ences of soil available nitrogen were significant between treat¬ 
ments in 2010-2011 (no data). Therefore, these results might be 
explained by nitrogen accumulated in soils during 5-years of 
annual fertilization, so that the available nitrogen was still higher 
than on control plots after nitrogen addition ceased. In a similar 
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study in a pine forest, soil NH 4 + gross mineralization rates on 
fertilized plots remained higher than on control plots, fourteen 
years after nitrogen addition ceased (Chen et al. 2006). 

A second explanation is that soil nitrogen availability was in 
accordance with the controls, but hysteresis or limitation of soil 
seed banks delayed the recovery of plant diversity. Shifts in re¬ 
source availabilities to match those of control plots might not be 
sufficient for recovery of the original community when hystere¬ 
sis affects communities characterized by resource-ratio competi¬ 
tive interactions (Suding et al. 2004). One study in a formerly 
fertilized field showed that short-tenn seed limitations were the 
dominant factors preventing recovery of plant diversity (Clark et 
al. 2010). 

Relative species richness, relative diversity and relative even¬ 
ness increased annually after cessation of fertilization. This 
demonstrates that the effects of fertilization declined and plant 
species were able to recover with each year after fertilization was 
stopped. However, because our study ended in 2011, we cannot 
estimate the time span required for species diversity to recover to 
levels recorded on control plots. Long-term experiments are 
needed to observe the dynamics of species diversity of sandy 
grasslands. A study with fertilization history similar to that of our 
study area estimated that grassland vegetation would need nine 
years to revert to the state of control plots after five years of 
nitrogen fertilization at 20 g N-m'^a' 1 (Mountford et al. 1996). 

Species composition changed on our study area after 3 years 
without addition of nitrogen. We attribute this to lower concen¬ 
trations of nitrogen in soils that altered interspecific competitive 
interactions after nitrogen addition ceased. Although the domi¬ 
nant species differed by treatment, the dominance of C. sativa 
and P. communis decreased from 59% in 2008 to 45% in 2011. 
This was probably because these species had higher nitrogen use 
efficiency under the nutrient enrichment conditions in 2008, but 
had reduced competitiveness when nitrogen availability declined 
after fertilization ceased. Community recovery is a complex 
process and is affected by many factors. Spiegelberger et al. 
(2006) studied mountain ecosystems and reported that resilience 
can be low in response to perturbations that substantially alter 
soil pH. Josef et al. (2009) hypothesized that the restoration of 
grasslands after cessation of fertilization was more rapid in pro¬ 
ductive grasslands than in unproductive grasslands. On our study 
area nitrogen is the main nutrient limiting net primary production 
(Yu et al. 2006). Inorganic nitrogen concentration in soils was 
significantly increased and soil pH was significantly decreased 
after five consecutive years of nitrogen addition (Zeng et al. 
2010; Li et al. 2011). We predict that the recovery of plant com¬ 
munity will require long periods of time in N-poor sandy grass¬ 
lands. 

The response of plant composition to cessation of nitrogen 
inputs differed by functional group. The important value of an¬ 
nual plants declined and treatments were similar in 2011. This 
was probably due to the lower dominance of annual-C 4 plants. 
Salsola collina and Digitaria sanguinalis disappeared in 2011. 
The annual-C 4 plants composition on fonnerly fertilized plots 
were consistent with levels on control plots in 2011. Biennial 
plants were the composits Artemisia scoparia and Artemisia 
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sieversiana. The differences of biennial plant importance were 
not significant between treatments in 2008-2011. But the impor¬ 
tance values of these species increased after nitrogen addition 
ceased. Compared to annual plants, perennial plants showed 
fewer signs of recovery. The importance values of perennial 
plants differed significantly between treatments during 
2009-2011. Two species that were present during the fertiliza¬ 
tion period, Lespedeza hedysaroides and Astragalus melilotoides, 
were not recorded in 2011, the third year after fertilization 
stopped. The recovery rate of species composition was slower 
than that of species richness, similar to the situation on prairie 
grassland, where relative species numbers, but not species com¬ 


position, recovered to control plot levels 13 years after nitrogen 
addition ceased (Clark et al. 2008). 

Vegetation height and coverage of fonnerly fertilized plots 
were greater than on control plots in 2011. Relative height and 
coverage declined over time. The recovery rates for coverage 
differed by functional group after nitrogen addition stopped. The 
relative coverage of annual plants was significant lower and the 
relative coverage of perennial plants was significant higher in 
2011 than in 2009. In contrast, the differences of relative cover¬ 
age of biennial plants between years were not significant (Fig. 3). 



Year 


Fig. 3: Changes of relative coverage among functional groups between treatments in 2009—2011. 


Conclusion 

Species richness recovered faster than species composition in 
Keerqin Sandy Lands after cessation of nitrogen inputs. Annual 
plants were more sensitive than perennial plants to the cessation 
of nitrogen inputs. We conclude that the recovery rate differed 
not only between community characteristics (species diversity 
and species composition), but also between functional groups 
within the same community characteristic. The different re¬ 
sponses were probably due to differences in N acquisition strate¬ 
gies and sensitivities to the cessation of nitrogen addition be¬ 
tween short-lived plants and long-lived perennial plants. 
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